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The Aharonov-Bohm effect is investigated in two-dimensional, single-terminal quantum rings in 
magnetic fields by using time-dependent density-functional theory. We find multiple transport 
loops leading to the oscillation periods of h/(en), where n is the number of loops. We show that 
the Aharonov-Bohm oscillations are relatively weakly affected by the electron-electron interactions, 
whereas the ring width has a strong effect on the characteristics of the oscillations. Our results 
propose that in those experimental semiconductor quantum-ring devices that show clear Aharonov- 
Bohm oscillations the electron current is dominated by a few states along narrow conduction chan- 
nels. 
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I. INTRODUCTION 

The Aharonov-Bohm (AB) effect is one of the most 
distinct physical phenomena which illustrates the impor- 
tance of the quantum mechanical phase. 1 ' 2 In the AB 
effect, a charged particle acquires a phase shift 



2ire/h / A • dr 



(1) 



from the vector potential A while traveling along the 
path 7. The existence of this phase shift can be verified, 
for example, by measuring the current through a quan- 
tum ring (QR) in a uniform perpendicular magnetic field. 
The electrons traveling along right and left arms of the 
QR gain a relative phase shift 
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ds = 27T$/$o, 



(2) 



where the surface integral yields the magnetic flux <I> 
through the QR, and $0 = h/eis the magnetic flux quan- 
tum. The resulting current is then a periodic function of 

The AB oscillation was first observed in 1960 by Cham- 
bers et a/., 3 and the first experiment with the period of 
<E>o in a QR configuration was made in 1985 by Webb et 
al. A who measured the magnetoresistance of submicron- 
diameter Au rings. Similar conductance oscillations were 
later found also in semiconducting QRs. 5 10 In addition 
to being an important concept in quantum mechanics, 
the AB effect has emerging applications in future com- 
puter technologies, random number generation, electron 
phase microscopy, and holography. 11 

In theory, AB oscillations in QRs have been studied 
both analytically and numerically. Analytic works 12 16 
have focused on one-dimensional (ID) rings where 
the electron path itself is not affected by the mag- 
netic field, and thus the situation is different from 
the experiments. 4 10 Numerical studies include (single- 
particle) wave-packet simulations in a ID system, 17 ' 18 



and very recently also in a two-dimensional (2D) 
geometry, 19 as well as 2D tight-binding calculations for 
conventional QRs 20 and graphene rings. 21 These numer- 
ical methods are able to take the effect of Lorentz force 
into account, in addition to the AB phase shift induced by 
the vector potential A. The largest effect of the Lorentz 
force is the asymmetric magnetic-field dependent arm in- 
jection of electrons which leads to the decreasing ampli- 
tude of the AB oscillations as the magnetic flux through 
the QR is increased. 17 ' 19 The asymmetric arm injection is 
qualitatively taken into an account in a recent S-matrix 
study by Vasilopoulous et al. 16 leading to a good agree- 
ment with numerical results in Refs. 17 and 19. In addi- 
tion to these efforts, the two-particle AB effect has been 
analyzed in rings where the combined path of indepen- 
dent electrons enclose the magnetic flux. 22 

The aim of this paper is to numerically study the 
transmission effects in a 2D semiconducting QR struc- 
ture in a static and uniform magnetic field. We use time- 
dependent density-functional theory which allows us to 
examine the real-time dynamics, the electron-electron in- 
teractions, and the role of the finite ring width (2D char- 
acter) in the same footing. We find that the width of the 
conduction channel is a critical parameter for the ampli- 
tude of the AB oscillations, whereas the electron-electron 
interactions have a relatively small effect on the oscilla- 
tion characteristics. The results give important guidance 
in assessing the transport dynamics in QR experiments. 

The rest of the paper is organized as follows. In Sec. 
II. we describe the QR structure of interest and the theo- 
retical model of the system. The computational methods 
used in this work are presented in Sec. III. The results 
for different number of electrons in the QR are presented, 
analyzed, and compared qualitatively with experiments 
in Sec. IV. The work is summarized in Sec. V. 



II. MODEL 

We consider a model for a semiconductor QR device 
fabricated in AlGaAs/GaAs structures. 5 10 Following the 
conventional approach in modeling the confined conduc- 
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FIG. 1: (Color online). External potential describing the 
quantum ring. The exact form of the potential can be found 
in Eq. (Al). Inset: Schematic image of the quantum ring 
potential, where the potential for the initial state (t — 0) is 
outlined by the dotted line. 

tion electrons in the material we apply the effective- 
mass approximation with the GaAs parameters, i.e., the 
effective mass m* = 0.067 mo and the dielectric con- 
stant e = 12.7 so- Hence, the energies, lengths, and 
times scale as E% = (m* / 'm ) / \e / 'so) 2 E h w llmeV, 
al = (e/eo)/(m* /m )a « 10 nm, and tg h/E% w 60 fs, 
respectively. In the results we use these effective atomic 
units unless stated otherwise. 

The time-dependent Hamiltonian describing our N- 
electron QR is given by 

N 1 

H = f + V ee + T> ext = J2 H fiV * + eA ( r 0f 
i=i 

e 2 N 1 * 
+ i^Sj^i+EWr*,*), (3) 

Z>J 1 J 1 1 = 1 

where the vector potential is given in the symmetric 
gauge, A = B/2{— y, x, 0), for the static and uniform 
magnetic field perpendicular to the QR plane. We ne- 
glect the spin effects by considering spin-compensated 
systems, so that the Zeeman and spin-orbit terms are 
omitted in the Hamiltonian. We construct a model for 
a single-terminal QR consisting of N = 1 ... 10 electrons 
and represented by the external potential v ext (r,t). As 
visualized in Fig. 1, the potential consists of the QR con- 
finement and the output terminal. Both the ring and 
the terminal have a Gaussian-shaped cross section and 
a tunable width (see Appendix). In addition, we have a 
linear potential describing a "bias" which is combined to 
a rectangular potential well at the end of the output ter- 
minal. A similar single-terminal geometry for a metallic 
QR (in ID) was considered already in Ref. (14). 



The static initial state is calculated for a quantum well 
located in the back part of the QR (see the inset of Fig. 1). 
Then, the QR and the terminal are suddenly opened up 
and the electrons are propagated (at t > 0) from the 
quantum well into the QR device. The opening is de- 
scribed by the Heaviside step function, which is in fact 
the only time-dependent part in the external potential 
(and in the Hamiltonian). Details of the potential can 
be found in the Appendix. 

Using a single-terminal device guarantees that there 
are no undesired effects of electron back-scattering at 
the input lead. To minimize the back-scattering at the 
output of the terminal, the length of the channel is rel- 
atively large (approximately three times the diameter of 
the QR). Moreover, we use a large potential well as a 
"sink" for electrons beyond the output of the terminal. 
The side length of the rectangular sink is larger than the 
length of the terminal. 



III. METHOD 

To examine the time-dependent many-electron sys- 
tem described by the Hamiltonian in Eq. (3) we apply 
time-dependent density- functional theory 23 (TDDFT). 
TDDFT is a practical, yet in principle exact approach to 
many-body dynamics, and it has received significant pop- 
ularity in, e.g., calculating electronic excitations. Quan- 
tum transport 24 can be considered as a recent field for 
TDDFT and it has been approached by both Green func- 
tion formalism 25 and finite-system time-propagation 26,27 
- the latter being the framework in this study. 

TDDFT rests on the Runge-Gross theorem 28 stating 
that for a given initial state there is one-to-one correspon- 
dence between time-dependent one-body density n(r,i) 
and the time-dependent one-body potential v ext (r,t). It 
means that a certain time evolution of the density is gen- 
erated by at most one time-dependent potential. Thus, it 
is possible to define an auxiliary Kohn-Sham (KS) system 
of noninteracting electrons moving in time-dependent ef- 
fective KS potential, such that the density of the nonin- 
teracting system is equal to the density of the real system. 
The KS potential consists of v ex t(r,t) (see above), the 
Hartree potential corresponding to the classical Coulomb 
interaction, and the exchange-correlation (xc) potential 
including all the indirect many-body effects. The xc po- 
tential is generally non-local in space and time and needs 
to be approximated in practice. Here we apply the com- 
mon adiabatic local-density approximation, so that we 
locally and instantaneously use the xc potential of the 
static and uniform 2D electron gas. We point out that 
the ring width is varied within a range that does not lead 
to the recently reported failure of the 2D local-density 
approximation in the quasi- ID limit of QRs. 29 

All TDDFT calculations are done with the octopus 
code package 30 published under the GPL license. The 
code has been built on the real-space grid discretization 
method which allows realistic modeling of QRs in 2D. 



3 



The calculations proceed as follows: the initial-state con- 
figuration is calculated by solving the KS equations for 
the initial external potential at t = and a constant, uni- 
form magnetic field B perpendicular to the plane. The 
eigenvalues are solved with the conjugate-gradient algo- 
rithm. After the initial calculation the KS orbitals are 
propagated with the external potential at t > repre- 
senting the full device while keeping the magnetic field 
constant. The propagation of the KS orbitals is done by 
applying enforced time-reversal symmetry approximation 
for the time-evolution operator U(t,to). After the time- 
dependent calculation the process is started again with 
a different magnetic field strength. We point out that 
the final propagation time is limited by undesired back- 
scattering effects. Therefore, to obtain a measure for the 
conductance, we monitor the electronic density up to a 
specific instant in time as explained below. 
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FIG. 2: (Color online). One-electron density in a QR with 
width parameter a — 0.5 (see Appendix) at different times. 
Magnetic flux through the QR is $ = 3$o- The symmetry of 
the density distribution is distorted due to the Lorentz force. 



IV. RESULTS 

As a measure for the conductivity in our finite system 
we consider the relative amount of transferred electrons 
at time t through the ring, 



JVr(*,t) = 1 



N ' 



(4) 



where TV = 1 ... 10 is the total (fixed) number of elec- 
trons, $ is the magnetic flux enclosed by the ring, and 
N(Q>,t) is the number of electrons in the QR at time t, 
which is obtained by integrating the electron density over 
a square which contains the QR. 



A. One-electron transport 

We begin with a single electron propagated through 
the QR. As an example of the dynamics of the system, 
the electron density with magnetic flux $ = 3<I>o is visu- 
alized at different times in Fig. 2. At t = the density 
starts to flow through the both arms of the QR. The 
densities in left and right arms collide near the output at 
t ~ 2 creating a standing wave that oscillates clockwise 
and counterclockwise in the system, while a part of the 
density is pumped out from the device. The asymmetric 
density distribution is caused by the Lorentz force. 

The transferred probability density at t = 1 and t = 2 
is plotted in Figs. 3(a) and (b), respectively. We have 
considered three different width parameters for the QR 
channels, a = 0.1, 0.5, and 1 (see Appendix), so that 
larger a corresponds to larger actual QR width. At t = 1, 
a smooth and regular AB oscillation with a period of 3>o 
is clearly visible in the conductance. The conduction 
minima and maxima occur at magnetic flux values cor- 
responding to phase shifts tt and 27r, respectively [see 
Eq. (2)]. An exception is the widest ring with a = 1, 
where we find a large, almost a half-period phase shift. It 
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FIG. 3: (Color online). Transferred probability density for 
a single-electron quantum ring at times (a) t = 1 and (b) 
t = 2 as a function of the magnetic flux $/$o- (c-d) Corre- 
sponding Fourier transforms of the situations in (a) and (b), 
respectively. Larger a corresponds to a wider ring. 



is due to transverse density variations and large scatter- 
ing effects in comparison with thinner rings. In addition, 
the amplitude of the AB oscillations strongly decreases 
as a function of a. This is due to the fact that the effect 
of destructive interference becomes smaller when there 
junction area increases. In other words, the left and right 
parts of the electron wave packet do not overlap as much 
as with a thinner output terminal. Increasing further 
the width of the ring and the output terminal eventually 
leads to the vanishing of the AB oscillation. Similar de- 
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structive effect due to the increasing QR width has been 
found by Pichugin et al. 20 

As shown in Fig. 3(b), at t = 2 the regular AB oscil- 
lation with the period of $o is distorted by oscillations 
with larger frequencies. The Fourier spectra of N r at 
t = 1 and t = 2 in Figs. 3(c) and (d), respectively, show 
that at first the major part of the electron transfer occurs 
at a period of $o in the flux frequency. At larger times 
the situation changes in an interesting way, namely, we 
find additional oscillations with AB periods of $o, ^o/2 
and $o/3. These oscillations arise from multiple trans- 
port loops in the ring: a part of the electron density does 
not flow out from the QR after traveling through one of 
the arms just once, but it continues to flow in the ring 
and can travel through an odd number of QR arms before 
exiting the ring. A similar phenomenon has been recently 
found in tight-binding 21 and wave-packet calculations for 
QRs. 19 Multiple transport loops have also been observed 
experimentally by Chang et al 10 

Next we analyze the multiple loops in more detail by 
considering the initial current split into two QR arms. 
The right hand part of the current receives a total phase 
shift 



4>R,tot&\(n r ) = ()> r + (n r - 1) {<t> r - <j>i) . 



(5) 



where n r is the number of times the current travels 
through the right arm, and cV (0/) is the phase shift 
obtained while traveling through right (left) arm once. 
Similar equation holds for the left hand part of the cur- 
rent 



^L.totalfaz) = <f>l + (ni ~ 1) {(j>i ~ </> r ) . 



(6) 



The total phase difference of the interfering currents can 
be obtained from Eqs. (2), (5) and (6) as 



2tt$ 



PR,total 



^L, total 



$o/(n r + m - 1)' 



(7) 



which leads to oscillation with periods of &o/n, where n 
is a positive integer. The amplitude of the flux frequency 
corresponding to the AB period of <po does not increase 
considerably at t = 1...2: the probability density trav- 
eled through only one arm is small at the output of the 
QR at t = 2. Proceeding further in time leads to more 
amplitude peaks for flux frequencies corresponding to pe- 
riods <&o/n. According to 2D wave-packet calculations by 
Chaves et al. 19 this effect is not visible in systems having 
smooth lead-ring connections. 

The effect of the Lorentz force is visible in the densi- 
ties plotted in Fig. 2. The perpendicular magnetic field 
favors the other arm of the QR and produces an uneven 
distribution of the density in the arms. The effect is not 
that visible in conductance graphs in Fig. 3 since our 
electrons start from rest. Szafran et al. report increasing 
maximum of the conductance oscillation together with 
the decreasing amplitude when the magnetic field is in- 
creased: the Lorentz force bends the path of the incom- 
ing electrons so that the probability of electron reflection 
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FIG. 4: (Color online), (a) Transferred relative number of 
electrons for a six-electron quantum ring at time t = 3 as 
a function of the magnetic flux <3>/<I>o- (b) Corresponding 
Fourier transform of the situation in (a). Larger a corresponds 
to a wider ring. 



at the input of the device decreases. The effect of the 
Lorentz force is much smaller in our simulations, since 
we do not model the scattering effects at the entrance of 
the QR device. 



B. Many-electron transport 

Next we consider the electron transport with several 
electrons. The relative transferred number of electrons 
for N = 6 and the corresponding Fourier spectrum are 
shown in Figs. 4(a) and 4(b). We find that increasing 
the number of electrons does not remove the regular AB 
oscillations, but the system seems to become more sen- 
sitive to changes in the width of the ring and the output 
terminal. The sensitivity appears as strong variations in 
N r as a function of a. They are due to complex electron- 
electron interaction effects, which, especially in the case 
of large a, have also significant transverse contributions. 
The thinnest ring with width parameter a = 0.1 shows 
signs of oscillations with phase difference of 7r emerging 
from the interference of scattered and unscattered elec- 
tron densities. In this case the scattering occurs predom- 
inantly along the transport channel. 

An explicit view on the effects of electron-electron in- 
teractions is given in Fig. 5, where we compare the non- 
interacting situation (solid lines) to the interacting one 
(dashed lines). The conductances are very similar. This 
is the case even if the radius of the ring is doubled [see 
Fig. 5(b)], so that the system is relatively more strongly 
interacting (note that the Coulomb energy scales as r _1 , 
whereas the kinetic energy scales as r -2 ). Although the 
enlargement of the QR notably increases the difference 
between the two situations, the qualitative picture re- 
mains; most importantly, the AB periodicity remains the 
same. However, the interactions seem to make the local 
minima at = n<&o l ess pronounced compared to system 
without interactions. This corresponds to a smaller am- 
plitude in oscillations with phase shift of it. Generally, 
the increase in N leads to the damping of the oscillation 
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FIG. 5: (Color online). Transferred relative number of elec- 
trons for different quantum rings with (dashes lines) and with- 
out (solid lines) interactions. The width parameter is a = 0.5 
in all cases, (a) N = 6, ro = 1. (b) N = 6, ro = 2. (c) 
N = 10, r = 2. 



amplitudes as shown in the comparison of the six- and 
ten-electron QRs in Figs. 5(b) and 5(c). 

Our results above show that the conductance is sensi- 
tive to both the width of the QR as well as the number 
of electrons confined in the ring. The effects are, how- 
ever, different in the sense that while increasing the width 
leads to drastic changes in the AB oscillations even for 
small N, the increase in N mainly damps the oscilla- 
tion amplitudes. In any case, it can be deduced that 
in semiconductor QRs containing dozens or hundreds of 
strongly scattering electrons the AB oscillations are un- 
likely to be visible. Hence, in view of the several exper- 
iments on semiconductor QRs where the AB oscillations 
are clear, 4 7 we may conclude that in a successful trans- 
port measurement the current in the semiconductor QR 
is induced by only the few highest states along a narrow, 
almost ballistic path. In metallic QRs (see, e.g., Ref. 4) 
the situation might be very different. 

To qualitatively assess the relative width of the elec- 
tron path compared with the QR radius we may consider, 
for example, the experiment of Fuhrer et al.J where the 
ratio between the width and the radius of the QR device 
is Ar/ro ~ 0.5. To obtain similar AB oscillation ampli- 
tudes as in the experiment for only a single electron, we 
need to set Ar/ro ~ 0.3 (here to estimate Ar we have cut 
the density profile at a few percent of its maximum) . This 
geometry already corresponds to a significantly thinner 
ring than the experimental one, and the difference be- 



comes even larger if TV is increased, since then the width 
in the calculation should be decreased to maintain the 
oscillation amplitudes. This qualitative analysis suggests 
that the actual electron path in the experiment is very 
narrow compared with the width of the device itself. 



SUMMARY 



We have investigated the Aharonov-Bohm effect in a 
many-electron two-dimensional quantum ring with the 
time-dependent density functional theory by modeling 
the discharging of a one- terminal device. We have found 
multiple transport loops leading to Aharonov-Bohm os- 
cillation periods of <&o/n = h/(en), where n is the num- 
ber of loops. The Aharonov-Bohm oscillations are rel- 
atively weakly affected by the electron-electron interac- 
tions, whereas the ring width has a strong effect on the 
characteristics of the oscillations. In general, the oscilla- 
tions are notably distorted when the conduction channels 
are made wider. Secondly, the increase in the number of 
electrons leads to the damping of the oscillation ampli- 
tudes. Our results suggest that in experiments on semi- 
conductor quantum rings where the Aharonov-Bohm ef- 
fect is observed, the actual electron path is narrow com- 
pared with the device size, and the transport is domi- 
nated by a few electrons close to the Fermi level. 



Appendix A: Quantum-ring potential 

The external potential in the Hamiltonian [Eq. (3)] is 
yiven by 



v ext (r,t) = 0(t)Vi(r) + (1 - 0(t))V o (r), 



(Al) 



where is the Heaviside step function, Vo is the ground- 
state potential (t = 0) and V\ is the potential used in the 
time propagation (t > 0). The potentials can be written 
as 



Vb(r) = (-Ve~\ = ) 

x0(x 2 +y 2 -(r o + d/2) 2 ) 

I _ e -(r /a) 2 V 

x6((rv+d/2) 2 -x 2 -y 2 )) 
x6(-x)6(y 2 c -y 2 ) 



<r /a) 2 
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and 

V 1 {v)=mm(-Ve {v/a)2 e{x),-Ve-^^ 1 ^^) 2 ^ 
x 6(L pipe + d/2 - x)6{x 2 +y 2 - (r + d/2) 2 ) 
^_fe-(^^) 2 -e-(-/«) 2 ) 

x 0((r o + d/2) 2 - x 2 - y 2 ) 

+ ^(-x/Lpipe - V/V'0(x - L pipe - d/2)), 

where V is the depth of the potential ring, V is the 
steepness of the linear potential, ro is the radius of the 
potential minima of the ring, y c is the cutoff length of 
initial quantum-well potential in the y-direction, a is the 
width parameter of the ring and the channels, L p i pe is 
the distance between the center of the ring and the edge 



of the potential well, and d is the grid spacing in the 
simulation box. 

Input parameters used in this work are V = 300, V = 
150, L p i pe = 10, ro = 1, ro = 2, a = 0.1, a = 0.5, a = 1.0 
and d = 0.005. In the propagation we have used time 
steps of At = 0.0002 . . . 0.0005. The length and width of 
the simulation box are 45 and 10, respectively. 
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